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Abstract 

Using very general and well established ideas of the statistical physics of macroscopic bodies, 
that is, of those composed of many degrees of freedom, we show how classical behavior of the 
center of mass motion arises from a fully quantum mechanical description of the dynamics of the 
whole body. We do not attempt to provide a rigorous proof of the latter statement, but rather, 
we show or, at least, indicate the hypotheses needed to obtain the purported result. Moreover, we 
neither attempt to deal with the “most general” physical situation and, instead, we concentrate 
on a stylized model of a small solid, yet macroscopic, that we shah call a “little stone”. The 
main hypothesis is that a macroscopic body can be decomposed into several smaller pieces, still 
macroscopic, that become statistically independent due to the short-range interaction nature of 
their constituent atoms. The ensuing main result is that the quantum distributions of extensive 
variables of the body become sharply-peaked. The center of mass variables are of this type and 
hence their dynamics is essentially classical. We point out the crucial role played by the external 
potential, in which the motion occurs, as the macroscopic agent that executes the “measurement” 
process of the center of mass. 
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INTRODUCTION 


Deriving classical behavior from a quantum one is a question that we all, practitioners 
of physics, have faced at some point of our careers. This is because we believe that the 
world is quantum and that classical behavior is an approximation. We are also aware that 
classical mechanics should apply to macroscopic bodies, such as balls, stones, houses, rockets, 
planets, et cetera. We know very well that at the atomic or molecular level, the dynamics 


must be described with quantum mechanics [l| , although there are instances in which, even 
at that level, classical mechanics works fine for the description of many physical properties. 
Therefore, there is always the question of when and how classical mechanics “emerges” from 
an otherwise quantum description. While there are related questions to the one posed here, 
that go under the names of “decoherence” and the “measurement” problem}^ |, our goal 
is not to dwell into concepts like those, but rather to simply make use of well established 
“facts” of the discipline of statistical physics j^, when explicitly applied to bodies composed 
of a very large number of atoms. True, we borrow those concepts, well defined in thermal 
equilibrium, and simply extend them to non-equilibium situations and, hence, we may leave 
dissatisfied more than one reader. But even if we cannot prove the main requirements, we 
believe that, at least, we point out what those requirements are, in order to obtain the 
classical limiting behavior of macroscopic bodies. 


As we shall elaborate in some detail below, the main hypothesis is that macroscopic 
bodies can be decomposed into smaller, but still macroscopic parts that can be considered 
as statistically independent. The physical reason behind this statement is the fact that atoms 
and molecules in macroscopic systems effectively interact via short-range potentials. As it 
can be shown[^, “short-range” means that the potentials decay as r“", with a > 3 and r 
the distance between the interacting atoms. Once one accepts this hypothesis, it follows 
that the (probability) distributions of extensive quantities are sharply-peaked, namely, that 
their fluctuations are I/'/N smaller than their mean values, with N the number of involved 
degrees of freedom. Hence, if iV ~ 10^° the values of extensive variables become accurately 
and uniquely determined, at any time. These statements are the basis of statistical physics, 
as clearly described in any textbook, such as Landau and Lifshitz Statistical 
What we will expose here is that the momentum and position of the center of mass of 
a (solid) macroscopic body are extensive variables and, therefore, their distributions are 
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sharply peaked. The obtained averages, or means, of those variables then obey classical 
mechanics, with “dissipative” terms that simply reflect the interaction of the center of mass 
with the very large nnmber of relative degrees of freedom of the body. The latter degrees 
of freedom do not need to obey classical mechanics; indeed, if the temperatnre of the “little 
stone” is low enongh, the relative degrees of freedom mnst be treated qnantnm mechanically. 
We insist that we shall focns on the model of a little solid stone becanse we want to “derive” 
Newton eqnations for its center of mass. If we dealt with a flnid, we wonld then obtain 
Navier-Stokes-like equations, a much technically harder endeavor, but that in any case could 
be realized with the same ideas here discussed; see also Ref. j^. An additional, and very 
important requirement to make possible the appearance of classical behavior, is the presence 
of an external anharmonic “macroscopic” potential. That is, for the different macroscopic 
subsystems of the little stone to be able to continuously “measure” the center of mass position 
and momentum, the relative and center of mass degrees of freedom must be coupled and, 
as we shall show, this can only be achieved if the external potential has anharmonic terms. 
Otherwise, the center of mass is always decoupled from the macroscopic number of relative 
coordinates, and one can no longer appeal to the statistical independence of its parts. It is 
interesting, therefore, to realize that the presence of macroscopic agents, namely, the external 
bodies that produce the external held, are ultimately responsible for the “measurement” 
and, in common but dangerous language, one may say that the external agent and the 
many relative degrees of freedom “collapse” the quantum evolution of the center of mass 
into a single, causal, classical evolution. This is not so unfamiliar. In equilibrium systems, 
the presence of a macroscopic conhning potential, such a box of rigid walls, is also partly 
responsible for the sharply-peaked values of their internal energy. 

Let us dehne a simple model of a little stone as a system of N identical atoms interacting 
pairwise in the presence of an external potential, 

N 


N ^ 

Pi 


^ = E ^ - ^il) + E ^ext(h;). 

i=i i<j j=i 


( 1 ) 


The “hats” on all variables is to insist that those are quantum operators. It is assumed that 
N ^ 1. Since we are assuming the temperature of the system is low enough to consider it as 
being solid, we should perhaps include more complicated interatomic potentials. However, 
the point is to make explicit that the interatomic potentials depend on relative coordinates, 
fi — fj. The interatomic potential is assumed to be short-range, namely that decays faster 
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than r~^ and their interaction range a is of atomic length scales, say a few Angstroms. The 
external held is “macroscopic” in the sense that its spatial variations are in length scales 
much larger than atomic ones; in other words, the quantization of a single atom in the 
potential Kxt(h) must give rise to energy levels that are separated in a scale much smaller 
than in any atomic energy contribution. It should be produced by other macroscopic bodies 
that do not feel the reaction of the stone, such as large coils, magnets, lasers or gravity. It is 
further assumed that all atoms feel the same external potential and, for simplicity, it should 
have a minimum. As it will be explicitly used, it must contain anharmonic contributions. 

There should not be any questioning into considering the system as “closed”, namely, 
that its only interaction with the “environment” is through the external potential. Under 
these conditions, the state of the system can always be considered as given by a density 
matrix p(t), at any time, that obeys Schrodinger equation, 

(2) 

provided a given initial condition /)(0). As it will be mentioned at the end, for systems 
with A^ 3> 1, in real life there is not much arbitrariness in the specihcation of the initial 
state. Therefore, although not essential, we shall assume that the initial state has a well 
dehned total energy, within the best experimental accuracy, that is, Et = (H) with variance 
= {H^) — such that SE Et- The average or expectation value of any physical, 
Hermitian variable of the system, A = A(ri,..., ... ,PAr), at any time t is given by, 

(i(()> = TV (p(()i) 

= IV (p(0)i(()) . (3) 

Because the interatomic interactions are short-range, the little solid stone can be thought 
to be composed of A/” macroscopic pieces, or subsystems, that interact with each other via 
the atoms at their borders only. The main assumption of this work is that those subsystems 
become statistically independent. This is of course true only approximately, yet this is the 
basic assumption of statistical physics. Let us denote by Ag, with s = 1, 2, 3,... , A/” variables 
that pertain to each subsystem. This means that Ag essentially depends only on degrees 
of freedom of the s-th subsystem. Using the full density matrix p, one can construct the 
distribution of values of any variable Ag by calculating all the moments of A at any given 
state, namely, calculating (A”) for all n = 1,2, ...cxd. This means that we can, in principle. 
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calculate W^As), the distribution of values of Ag in the given state. Statistical independence 
means that the distribution of all variables Ag, with s = 1, 2, 3,..., A/" is, 

hh(Ai, ^2,..., Am) ^ >V(Ai)>V(A2) ... W{Am). (4) 

While this is true in equilibrium, we extend it to “arbitrary” states, realizable in typical 
situations. 

Of essential importance are the variables which are extensive or additive. That is, those 
that obey, 

Ai-m = Ai + A 2 + ... + Am (5) 

where Ai_m is the value of the same variable, but for all the subsystems 1 to M together. 
These type of variables commute with each other. The main property of extensive variables 
is that they scale with the number of atoms of the corresponding subsystem j^. In particular, 
the average value of Ag is proportional to Ng, 

{As) ~ Ng, ( 6 ) 

and, therefore, (Ai_m) ~ W + A ^2 + • • • + Nm. The key result is, now, that the variance of 
an extensive variable is also extensive. Let us see. Consider the variance of Ai_m, 

MM MM 

(ip„) - = (EE-4T.) - 

j=l k=l j=l k=l 

M M 

= E(T> + - E{h>" - E(h)(4)- (7) 

j=l j^k j=l j^k 

We use now the fact that the subsystems are statistically independent to yield, 

{AjAk) = {Aj){Ak) for j ^ k. (8) 

Hence, we obtain, 

M 

{AIM - = E ((h') - ■ (9) 

i=i 

Namely, the variance of an extensive variable is also extensive and, therefore, it scales with 
the number of atoms in its subsystem, 

{AD - {Ag)^ ~ Ng. (10) 

As a matter of fact, all of the cumulants ((A^ — (As))"'), with n = 2, 3,4,..., are all extensive 
and, thus, all scale with Ng, i.e. ((A^ — (A^))") ~ Ng. 
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Since the fluctuation, or width of the distribution W(As), is proportional to the square- 


1/2 




, then this scales as SAg ~ On the other 


root of the variance 6As 


hand, as the mean value is at or very near the peak of the distribution 'W{As), then, 
the result 5Aa/{As) ~ l/y/Wl implies that W{As) is sharply peaked at the mean value. 


if S> 1. In other words, Ag essentially does not fluctuate and, effectively, looses its 


“quantum” property of “being” a superposition of more than one state: at every time, it is 
just in “one state”, its mean. This is not quite true, the system can still be in many different 
states, however, the density of states is so dense that, up to a precision 1/y/N^, one cannot 
experimentally distinguish among them. 

One may believe that there are many sensible physical variables that are extensive. This 
is not the case. As a matter of fact, the only extensive mechanical variables that can 
be constants of motion in a closed system are the total energy, momentum and angular 
momentum!^. In a typical system studied from a thermodynamic point of view, both the 
angular and linear momenta can be taken as zero since the system is assumed to be at 
rest in the lab. These conditions give rise to the preponderance of the internal energy 
in determining the states of systems in thermal equilibrium. Moreover, if the statistically 
independence of macroscopic subsystems holds true, this right away indicates why one can 
always associate well defined internal energies to macroscopic systems even if not in thermal 
equilibrium. One may say that the different subsystems “measure” among themselves their 
internal energy at any time - barring exceptional states that one may prepare in the lab. At 
the same time, this shows why it is difficult to “prepare” an otherwise closed system in a 
state that does not have its energy sharply peaked. With the previous ideas, the classical 
behavior of the motion of the center of mass follows very simply, as we now exhibit it. 

The degrees of freedom of the little stone, denoted by (fq,..., r7v;Pi, • • • ,'Pn)j can be 
separated into center of mass and relative coordinates. These are. 


- fj+i J = 2,..., A^ - 1. 



( 11 ) 


The corresponding conjugated momenta are. 


N 


p = T.Pi 


2 = 1 
N-1 


ft- = y Kjtpk j = 1,2,..., iv -1, 


( 12 ) 


k=l 
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where the symmetric {N — 1) x (TV — 1) matrix Kjk is equal to, 


^jk ^^j,k ^j,k—l 


(13) 


With this, we can write the Hamiltonian as, 

p N / ^ N—1 \ 

» = ^ + + (M) 

i=l \ k=l / 

where Hy-ei is the Hamiltonian of the relative coordinates (rfc,p^), k = 1, 2,..., — 1; it 

looks cnmbersome and so we do not write it. The qnantities are the coefficients of the 
inverse of the transformation of coordinates given by ffTT]) . not explicitly written, bnt with 
the next two nsefnl properties: 

N 

Eaf = 0 V j = l,2,...,iV-l 

i=l 

N 

^afaf = Kp. (15) 

The relevant point we want to make is that the relative degrees of freedom are deconpled 
from the center of mass coordinates {R,P), except for the external potential term, last in 
Eq. ffTT|) . It is this last term the responsible for the “measurement” of the center of mass 
position and momentum, within a very narrow latitude. 

Consider Erst two very common external potentials, one, the held produced by gravity, 
and other, a harmonic oscillator. 


= mgz 




(16) 


with z the coordinate in the direction of gravity. Note that m is the atomic mass. Using 
this in the last term of Eq. ffTTl) yields, 

N / ^ N-1 \ 


E Kit ^ + E 4* K = NmgZ 


i=l 

N 


k=l 

N-1 


N-1 N-1 


E Kxt° ( ^ + E 4 k ) = ^Nmu^R + -mu'^ E E ' K, (17) 

/ j=l k=l 


i=l 


k=l 


with Z the coordinate of R in the direction of gravity; to derive these expressions we used the 
identities ([15]). Note that in both cases the relative and center of mass motion are completely 
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decoupled. It is also of interest to highlight the known result that the mass associated to 
the center of mass coordinates is the total mass of the stone, Nrrt] it is of interest because 
it yields a “macroscopic” potential for that coordinate. It is hard to imagine a macroscopic 
stone oscillating at frequencies u of the order of megahertz or higher. So, we must consider, 
realistically, much smaller frequencies. This is in accord with the statement that the external 
potential must be macroscopic. From Eqs. flTTI) we hnd that, in order to couple the center 
of mass and relative degrees of freedom, the potential must have a n h armonic terms. In the 
case of a stone falling in the held of gravity, the presence of the hoor of the lab will do it: 
while the stone is falling, the relative motion does not feel gravity. However, during the 
collision of the stone with the hoor, the relative degrees of freedom couple to the center 
of mass motion and the kinetic energy of the center of mass is transferred to the relative 
motion, heating the stone up. 

Now, considering a generic anharmonic potential Kxt(fO) with an absolute minimum 
to make it simpler, and with spatial variations in macroscopic scales (i.e. with energy 
transitions much smaller than atomic ones), we can make an expansion around the center 
of mass coordinate i?, 

N / ^ N-1 \ 'I 

E pi + E 4“’?-= “ Nv,,t(R) + (18) 

i=l \ k=l J ^ OKaOKp 

where a and [3 are indices over the three spatial coordinates and summation is implied. The 
tensor operates on the relative variables only. It is given by, 

N-l 

h? = E KitijAn (19) 

fc=l 


and it is related to the moment of inertia tensor of the body. The hrst term in 0181) is the 
external potential felt by the center of mass only and the second one is the coupling between 
relative and center of mass motion; clearly, 14 xt(-R) must be anharmonic, otherwise there 
is no coupling. There may be situations, the collision of the stone with a floor mentioned 
previously as an example, where the above expansion may not be possible to perform. We 
shall assume the expansion is allowed just for purposes of exposition. 

Let us now write Heisenberg equations of motion for the center of mass, j = x,y, z, 


dt ^ 

±p 
dt ^ 


Nm 


= -N- 


d 


dR. 


■Kxt(i?) - ^x: 


1, 

2 ' 


(rel) 

'a/3 


d 




dR^ OR^dR^ 


-Kxt(i?). 


( 20 ) 
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Similar equation can be written for the relative variables and those are coupled to the above 
ones through the last term in Eq. flT8|) . 

Once again, for linear or quadratic external potentials, the last term in the second equa¬ 
tion of fl2U]) vanishes and the motion of the center of mass becomes decoupled from the 
relative coordinates. In such a situation the density matrix of the system can be written 
exactly as, 

P{t) = pCM{t)pve\{t), (21) 


where pcM and prei are the center of mass and relative degrees of freedom density matrices 
respectively. In this case, although the motion of the center of mass corresponds to a 
particle of macroscopic mass Nm, its dynamics is fully quantum and no classical limit can 
be obtained. 

However, if the center of mass and relative degrees of freedom are coupled, by the last 
term in the second equation of fl20|) . then. Erst of all, the full density matrix does not 
factorize as in Eq. fl2l]) . i.e., p{t) ^ /3cM(t)prei(t). However, it is this coupling what allows to 
obtain the classical limit of the center of mass motion, as we now argue. We recall, first, 
that both the center of mass position and momentum operators can be expressed in terms of 
the original positions and momenta of the iV particles of the sytem i = 1,2,..., N, 

see Eqs. (mi and (|T^ . The inverse transformation indicates that the original variables 
are functions of center of mass and relative variables. But since the latter are coupled, 
then, all variables become dependent among them. Now it enters the particularity of the 
model at hand, which is supposed to represent a little, solid stone. That is, because of 
the assumed short-range interactions, the original variables can depend on all of them only 
through collisions with their neighbors in the solid. This means that we can consider that 
the stone can be separated in even smaller pieces, yet macroscopic, such that the interactions 
between neighboring pieces occur only at their borders. Under these conditions we assume, 
once more, that these little pieces become statistically independent. 

To be a bit more precise, we can write the center of mass variables as. 


R 


1 

Mm 


N 




p = Y'Ps 


S=1 


( 22 ) 


where (TZstRs), s = 1,2,... ,M, are the center of mass variables of the s-th piece. Mlg are 
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the masses of the pieces. Since the pieces are macroscopic, the number while arbitrary, 
is not macroscopic and depends on the size of the stone. 

Due to all previous discussions, the variables of the centers of mass of the pieces are 
“weakly coupled” and, therefore, typically, are statistically independent. At the end of the 
paper we discuss the meaning of “typically”. Under this result (or assumption) and because 
of the relations given by Eq. fl22|) . the center of mass variables (i?, P) become extensive. As 
we have already shown, this implies that their fluctuations 5Ra and 5Pa scale as A^^A 
therefore, the (measurable) distribution of values of center of mass position and momentum 
are sharply peaked. 

Let us dehne the mean values of the center of mass variables as, 

m = {kt)) 

m = ikt)} ( 23 ) 


where the expectation value is taken with the full density matrix of the system p{t). Then, 
the equations of motion for the mean values are, 

P. 


R^= 


d 
d 


Nm 


dt 


P-, = -ivArrex.(fl) - i 


dR^ 


d 




dR^ dRr,dR, 


-Ve^tiR). 


(24) 


L'-y ^ KyJ. Liy KyJ~h^ 

In writing this equation we have assumed, for simplicity, that at this level of approximation 
the center of mass and relative coordinates are also statistical independent. Although we 
expect this condition to be true in typical states, it is not really necessary to impose it 
but certainly makes the argument much easier to follow. However, what we claim is truly 
correct is that R is indeed extensive, in the sense described above, in order to obtain the 
very important property, 

(Ue.t(^)) ~ Uext({^)). (25) 


Equation fl24)) is the main result of this note. It indicates that, indeed, the means of center 
of mass position and momentum obey Newton classical equations of motion. But since we 
can drop the adjective “mean” because the fluctuation and higher cumulants are negligible, 
we simply say that the center of mass variables do obey Newton equations. While Eqs. (I2T)) 
look like “simple” Ehrenfest equations of motion, we stress that their validity necessarily 
and explicitly require that the body is macroscopic. This fact is reflected in the last term 
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of the second equation in (|2^ . This term, as we analyze below, is of a “dissipative” nature 
in the sense that it is responsible for the transfer of energy between the center of mass and 
the relative degrees of freedom. Even if the statistical independence of center of mass and 
relative variables were not accurate, their coupling would result in a dissipative term for the 
center of mass. 

For purposes of exposition we have not insisted in the time evolution of the relative degrees 
of freedom. It is clear that we can make further assumptions for their study. For instance, 
one can assume that the atoms in the solid vibrate around their equilibrium positions. This 
would yield a further separation of relative variables into “rigid body” rotations and relative 
small oscillations that, in general, would be coupled. In the appropriate limit, the rigid body 
degrees of freedom should also obey classical equations of motion, coupled to the center of 
mass and the relative vibrations. We believe that the main point of this note does not need 
to make this further analysis explicit. 

To conclude the argument, we can now study the time evolution of the energy of the 
center of mass, dehned by the operator. 


Ecu — 


P 


2Nm 


+ ArKxt(^). 


(26) 


By itself, it is the Hamiltonian of an “isolated” particle of mass Nm in the presence of the 
external potential NVext- It is of interest to note that while it is fully quantum, it refers to 
a single macroscopic “particle”. The time evolution of Ecu is given by. 


d i V 

wt-E'cm = w H-,Ecm 
at h L 


(27) 


with H the full Hamiltonian, given by Eq. Using the assumed expansion of the external 
potential given by Eq. flT8|) . one hnds, to the same order of approximation. 


1 

dt 


E, 


CM 


j-(i-el) 

2Nm 


d 




dR^dR^dR 




fi 


(28) 


Dehning Ecuit) = Trp(t)EcM, and using the statistical independence of the pieces of the 
stone, we obtain 


Wt 


E, 


CM 


1 

2 



d 

dR^ dRadRjs 


VextiR) 


dR^ 
dt ’ 


(29) 


namely, a clear dissipative behavior that leads to an eventual mechanical equilibrium state 
where the energy of the center of mass of the stone will no longer change. It does not directly 
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imply that the energy of the stone is dropping all the time from the initial state. There 
may exist intervals of time of decrease or of increase of energy, depending on the details 
of the initial condition and the potential, but given that there are a very large number of 
relative degrees of freedom, iV 3> 1, we expect from general considerations of the irreversible 
behavior of macroscopic bodies, that eventually the energy of the stone will tend to a 
minimum. In this situation, all the extra initial energy of the stone will be released to the 
relative degrees of freedom. The stone, in turn, will end up with a higher temperature than 
initially. It is also interesting to observe that when we do classical mechanics theory for the 
motion of a stone, we usually neglect the dissipative term in Eqs. (^^ and fl25]) . thereby 
assuming that the motion of the center of mass in the external potential is conservative. 
While we do know that this is never true in real life, we hnd from the present argument, a 
glimpse of how the dissipation is always there and, furthermore, that its origin can be traced 
back to quantum mechanics. 

The present “derivation” or exposition of requirements, for the attainment of classical 
motion of the center of mass of a macroscopic stone, has certainly hidden the problems 
of “decoherence” and “measurement”, typical of all discussions of QM when applied to 
macroscopic bodies}^. To begin, we have avoided to discuss the processes, and the time they 
take, that lead an arbitrary initial state to one in which statistical independence has been 
reached; this is the problem of decoherence. However, we claim, “typical” initial conditions 
in the lab, for a real stone of atoms are already in states where statistical independence 
is essentially valid. One would have to work very hard to devise a procedure to prepare 
a stone in an initial state, say, in which its energy was not already sharply peaked. It is 
important to observe that the part of the system composed of the relative degrees of freedom 
can be at very low temperatures, such that a quantum many-body treatment would still be 
required. That is, the statement of the emergence of classical behavior refers to the center 
of mass degrees of freedom only. Yet, we do admit, we have not shown how decoherence 
actually occurs. The second issue, that of measurement, we have already alluded to it. That 
is, once we accept that statistical independence of the macroscopic little pieces that conform 
the stone sets in, these pieces “measure” among themselves all their extensive variables, 
giving them sharp values at all times. And we recall once more, the crucial role played by 
the macroscopic external potential in determining the extensive statistical property of the 
center of mass variables. But again, this nothing but the essence of the success of the theory 
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of statistical physics. 
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